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ON THE CONVERGENCE OF SPACE-TIME DISCONTINUOUS 
GALERKIN SCHEMES FOR SCALAR CONSERVATION LAWS 

GEORG MAY* AND MOHAMMAD ZAKERZADEH* 


Abstract. We prove convergence of a class of space-time discontinuous Galerkin schemes for 
scalar hyperbolic conservation laws. Convergence to the unique entropy solution is shown for all 
orders of polynomial approximation, provided strictly monotone flux functions and a suitable shock¬ 
capturing operator are used. The main improvement, compared to previously published results of 
similar scope, is that no streamline-diffusion stabilization is used. This is the way discontinuous 
Galerkin schemes were originally proposed, and are most often used in practice. 
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1. Introduction. In the present paper we consider space-time discontinuous 
Galerkin (DG) schemes for the scalar conservation law 

d 

(1.1) ut + '^Mu)a;i =0, in (0, oo) X =: 

(1.2) u(0, •) = Uq, in 

where / = (/i,..., /d) is a smooth flux function. We assume that uq € L°“(R'*) with 
compact support. A classic result states the existence of a unique entropy solution in 
L°°(R(i^^), still with compact support mm Recall that a weak solution to dm is 
an entropy solution if it satisfies, in the distributional sense, 


d 

(1-3) r]{u)t + ^ Fi{u)xi < 0 

i^l 

for any convex entropy function ri{u) and associated entropy flux function F = 
(Fi,..., Fd) that satisfy the compatibility condition F/ = f^ri' for i = 1,..., d. 

We prove convergence of a class of DG schemes to the unique entropy solution 
of (11.11) , (11.21) . Using DiPerna’s theory of measure valued solutions [7] , we can assert 
strong convergence in for 1 < p < oo, provided that the approximate solution is 

i) uniformly bounded in F°“(R(J+^), 

ii) weakly consistent with all entropy inequalities, 

iii) strongly consistent with the initial data. 

These statements will be made more precise below. It is worth noting that, for Cauchy 
problems of the type (EH), (EH), one may use Szepessy’s extension |19| to DiPerna’s 
theory, to replace the first criterion with the weaker requirement that the approximate 
solution be bounded in F°“(R+; F^’(R‘^)). We comment on this distinction below. 

Convergence results along these lines have been obtained previously for both 
continuous Galerkin (CG) and DG finite-element methods [TH [T31 [TO] ■ The schemes 
that were considered in the above references, both CG and DG, have in common 
that they use streamline-diffusion stabilization in conjunction with a shock-capturing 
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operator. At least in the case of DG methods, streamline-diffusion stabilization is 
not commonly found in practical implementations miiiiiin]. For linear hyperbolic 
problems, the analysis of DG schemes has long ago made the transition from such 
schemes that use streamline-diffusion stabilization m to schemes that do not [ 12 ]. 
However, the techniques used for linear problems have no direct extension to the 
nonlinear case. To the best of our knowledge, a proof of convergence, similar in 
scope to that presented here, is not available for DG schemes applied to nonlinear 
conservation laws, unless streamline-diffusion stabilization is used. (In addition to 
the references cited above, we also mention the error estimates provided in |^.) 

A proof of convergence exists for GG schemes without streamline-diffusion stabi¬ 
lization, using polynomials of degree q = 1 m- In the present paper, convergence 
of a space-time DG scheme for ini), (na), similar to that of m, is proved for 
uq G with compact support, and any degree <? > I of polynomial approxima¬ 

tion. The scheme does not use streamline-diffusion stabilization, but only a suitable 
shock-capturing operator. For our proof, we use DiPerna’s framework [7], requiring 
a uniform L°°-bound on the approximate solution, as opposed to the modification 
due to Szepessy m- The reason is that, firstly, showing the existence of such a uni¬ 
form bound is an interesting result in its own right. Secondly, the existence of such 
a bound facilitates the proofs in this paper, and thirdly, it allows the extension to 
intial-boundary value problems. (See, for example, the results in [20l [Tn] .l 

The paper is organized as follows: In section [2] we introduce the space-time DG 
scheme, and define some notation. In section [31 using a suitable test function, we 
re-write the scheme in a convenient form from which all the remaining results can 
be derived. We prove a uniform bound of the numerical solution in L°“, and the 
weak entropy consistency of our DG scheme in sections Sj and jS] respectively. Noting 
consistency with initial data, we are then in a position to conclude that the DG scheme 
converges to the unique entropy solution. 

2. Space-Time DG Scheme. We define the space-time coordinates x = 
(xo, xi, ..., Xd) with xo = t. It will sometimes be convenient to write x = (ii, ..., Xd) 
for the spatial variables. For arbitrary T > 0, consider a sequence of time instances 
0 = to < ^1 < ■ • ■ < tN = T, and define corresponding time intervals = (t„,t„+i). 
Let Td = {«} be a subdivision of the space-time slab /„ x into disjoint simplices. 
Define % ■= [jn^oV,"’ and let h := sup^g^-^ / ik , where := diam{K). (We may 
assume h < oo.) Convergence analysis will be carried out for shape-regular families of 
triangulations {Th}, which means there is a constant C > 0, independent of /i, such 
that 


( 2 . 1 ) 


1 / h^\dK\ ^ ^ 


Vk e Th- 


We define the approximation spaces 


N-l 

(2.2) K), = {u : vU G G m , Vh,^ = H 

n—0 

where is the space of d -I- 1-variate polynomials of total degree q. In the following, 
we assume that q > 1. (The case g = 0 is dealt with in Ref. |4].) Owing to compact 
support of the initial data, and the hnite domain of dependence, we can assume that 
u = 0 for large \x\. 
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The scheme which we analyze is defined by determining Uh G Vh^q^ such that 

(2.3) Af{uh,v) =M^^{uh,v)+Af^'^{uh,v) =0, \/v &Vh,q. 

Introducing the shorthand notation ^ := SugT"’ define the semi linear 

form 

(2.4) V ■ f{uh)vdx + J (^f - v ds^ , 


where / = {u,f{u)) is the space-time flux. On dn we define the interior trace values 
as = linij^Q- Uh(x + eiy), where ly = {vt, ■ ■ •, is the outward pointing 

normal on dn at x. We denote as Ke G 7h the element separated from k on e C 
dn, and introduce the numerical flux function / = defined to be 

Lipschitz-continuous, strictly monotone, and consistent with /. On interfaces where 
ly = (±1,0,..., 0)^ we assume that the numerical flux reduces to pure upwinding (i.e. 
in time). In this case, for any 0 < n < — I, we can re-write the boundary integral 

(2.5) 

^ f (f-f{ui^'>)-iy)vds+f (u)('-u;;')u+dx. 


where := {t„} x and 9 k* is the ’’internal” boundary of k G 7)^", consisting of 
points that are shared by at least two elements in x R'’*. Furthermore for x G R((, 
we define v'^ := lime_>ou(t„ ± e,x). 

We augment the semi linear form (12.4p with a shock-capturing operator 


( 2 . 6 ) 


J\f^^{Uh:V) : = 




e{uh)Vuh ■ Vvdx, 


with the viscosity parameter 


(2.7) e{uh)U := C, 


LN-f(.Uh)\dx + Jg^\u^ 


(«) 


ds 


f dx 

J K 


G Th, 


where ^ < /3 < 2, and Cg is a positive constant of order unity. (In the following we 
set Cg = I.) 

We define the norm 


and assume that ||/'||oo < OO0 This implies that 

(2.9) |V./(u)| < ||/'|U|Vu|, VvGVh,q. 


(2.8) ||/'||oo:=sup|/'(0|, |/'|:= 


In the following C shall denote a generic positive constant, independent of Uh and /i, 
but not necessarily the same at each occurrence. 

^This is not a severe restriction, because the true solution to problem l ll.lll . 1 11.21) satisfies a maxi¬ 
mum principle. In order to comply with the growth restriction, we may apply a smooth modification 
to f(v) for V outside the interval [min 2 ,^j[|d uojmax^^jgd uq]. 
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3. The Canonical Entropy Test. All results in this paper can be obtained 
by considering in (12.31) test functions of the form v = 4>). Here, tt/j is a suitable 

projection operator into the space Vh,q, is a smooth non-negative function, and 77 is 
a convex entropy. Following na, we define iTh using an iJ^-projection such that for 
all K G 71 


(3.1) 


j V{Trh{r]'(l))) ■ Vvdx = j V {r]'(j)) ■ V v dx, \/vG'P'^{k) 

/ •nh{'n'4>)dx = / r]'(j)dx. 

J K J K 


Using V = Trh{r]'4>) in (12.31) yields 


(3.2) 


JV^'"{uh,r]'(j)) +Af^^{uh,TTh{v'(l>) - +N'^^{uh,v'(l}) = 0 . 


Note that Af^^{uh, TTh{r]'(j))) = which is a consequence of the definition 

of the iJ^-projection (13.11) . We define the space-time entropy flux F := ( 77 , Fi,..., Fd). 
Then F satisfies the extended compatibility relation 

(3.3) i?/ = AV, 7 = o,...,d. 

The first term in (13.21) may be rewritten as 

■A/''°'^(uh,7V) = XI (/ ■ F{uh)<l)dx + (^f - r]'{u[^'’)(j)ds^ 

N-l 

+ J2 i^h -KW{u^)(pdx, 

where the boundary integrals are split as in (12.51) . Integration by parts yields 

(3.4) = - [ F{uh) ■'V(j)dx ■ ^(pds 

JR^ n.K 

+ X / / Tj{uj;)(j)dx 

n,K •'9 k, •'Rfr 

- ? 7 (uo)'('dx-h V / {r]iuj;)-rj{u+) + {u^-uj;)r]'{u^))(l)dx, 

“'K'" n=0 

where := ( 0 ,T) x has been defined, and we have equated = Ug for n = 0 . 
(In practice one might use a suitable projection for the initial conditions, but this is 
of minor importance here.) Using (13.31) . we have 


X/ ■I'c^ds 


1 


' Ok* 


/ 


{F(n'F)-F{uF'))-yt>ds 




2 J Ok, 


,(«) 

‘'h 




= X 


1 


f ■I'T]' 


> Ok,* 






„(-) 


,Ge) 


/(O ■ vri"{Cl dC \ <?^ds. 


(3.5) 
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After rewriting the second surface integral over i9k* in (13.4p to include contributions 
from both sides of the interface, and noting 

(«) 

we can substitute expression (ixa into (EH) to obtain 

(3.6) N^^{uh,ri'(i)) = - f F{uh) -Vcjidx + ( ri{u'^)(j)dx - f r]{uo)(j)dx 

jR'i 

N-1 „ 

+ Z! / {v(^h)-v{ut) + {u+-u-;;)r]'{u+))cl)dx 

_n 


n—0 ' 




' 5 k* J u 


.Oe) 


(«) 


(/(O - v- pj ?7"(0<(>d^ds. 


For later use, we note that, for any convex entropy, and non-negative function tp, 

(3-7) {viu^) - viul) + {u-l^ - uj;)r]'{u-l^)) </> > 0. 

Furthermore, by assumption, / is strictly monotone, which implies the E-flux property 
in the sense of Osher m, and so 


(3.8) 


L («« 


^ r]''{^)<j)d^ > 0. 


4. Stability. First, we state a discrete L^-stability result. Use r] = and 
(j) = 1 in EH, and observe that for this choice r]'{uh) € Vh^q- Consequently the 
second term in (13.21) vanishes. From (12.61) and (13.6F we obtain 


„(-e) 


E 


e{uh)\S/uh\^ dx+i 1 I 

: ^ JdK,^ 

N-1 


(/(C ds 


(4.1) 


+ EJ/ + I f K)^d2; = i/' uldx. 

n=0 ^ ^ ^ JRd. 


Therefore, noting (13.8L clearly we have discrete L^-stability in the sense 

(4.2) ||M/i(-,i)v)IU=(R‘^) ^ ll^o||L2(R,i), -(V=l,2,- 

At the interfaces R)(, equation (14.11) explicitly states an L^-bound on the jumps of the 
solution. For a strictly monotone numerical flux, a bound on the interface jumps of 
the solution Uh at is also implied [5], i.e. there exists a constant C > 0 such that 

r ^ r 

(4.3) / if{0-v-f{ul^\ul^‘'>;v))d^ds>C ds. 

J dK.» J dK* 

Now take rj = with an even integer p > 4, and (^ = 1 in (13.21) . Using (13.6F 
and the inequalities (13.7L (13.81) . we obtain immediately 

(4.4) 
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The term involving the projection error can be written 

f V ■f{uh)iTrhiul~^)-ul~^)dx 

n,K 

JdK J 

The i7^-projection (13.11) is a simple Neumann problem. For a smooth enough function 
V, we have the standard error bound ||u — '!Thv\\L°°(ui) < Ch‘^^^\\v\\\Yi+'i-,°°{K.) for lo = k 
or w = dK. Setting v = where w € Vh,q, we can further rewrite the projection 

error (cf. [20] ') 


(4.5) 






2^9+1 IU„l|P-3 


Using Lipschitz-continuity of the numerical flux, and the definition of the artificial 
viscosity (l2Jt . we thus obtain 

n,K 

Split the triangulation into Tf^ := {k G Th ■■ \\uh\\L'^(K) > 1} and ■= {k G % : 
\\uh\\L’=°(K) < 1} to obtain 

\Ar^^{uh,Trhiul~^) -ul~^)\ < 

f ^K)ll«^llL~%)ll^«'‘llL“(«)da;+ f 4uh)\Vuhfdx 

where the standard estimate 


^6r< • 


(4.6) |N|||Vu„||i,.(,) < ||Vu„||i.(,) 

has been used. To proceed, we need the following lemma: 

Lemma 4.1. For a shape regular triangulation Th = {n}, there is a h-independent 
eonstant C{q,p) sueh that \/k gTh 

(4.7) 

[ lblli,~%)l^^l^dx < (7(9,^) f Vv ■V{TTh{v^~^))dx, V G'P‘’{k),p = 2,4,--- . 

J K J K 


The proof is given in the appendix. This lemma was proved in [18] for linear elements. 
In (TU] the authors give a proof for v G V^^k) and a Lagrangian interpolation 
operator. In the appendix we present a novel proof for the lemma, which is simpler 
and applicable to the iJ^-projection (13.11) . Using Lemma [4.11 L^-stability dlTD, and 
(gH), there follows 


|A7^^(u,„7r,«-') - <-')! < 
(4.8) 


C'(<7,p)/r2-V+'E 


£{uh)VUh-V{'Kh{vfh ^))dx. 


Finally, substituting this result into (14.4L we have 
WMtN^ ■)\\lp(Rd)+p{^-C{q,p)h^-^p'^+^)M^'^{uh,TTh{ul~^)) < 
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Clearly, in order to keep the norm bounded or all n = 1,..., iV, it is required 
that 


< oo, 

P-C((7,p)/i2-V+^ > 0. 


From these conditions we get 


(4.9) 


^ < min{ 




}, 


and note that as h —> 0 this will still permit us to let p —oo eventually. 

It is worth pointing out that so far the presence of streamline-diffusion stabiliza¬ 
tion (or lack thereof) is clearly not essential. In previous proofs [TSJ [201HH] > the control 
over the residual provided by the streamline-diffusion stabilization becomes significant 
when extending the discrete stability to a uniform bound. We now establish a uniform 
bound in without using streamline-diffusion stabilization: 

Proposition 4.1. Let Uh be the solution produced by scheme (1231) . There exists 
a constant C > 0 such that 


(4.10) 


ll^/t||L“([0.T]xR'^) ^ ^ 


Proof. Fix t G [t„, t„+i] for some n = 0,..., — 1. For k £ Tff we define C n 

as the subset of k with first coordinate xq = t G [tn j t] , and analogously we define 
dnf C 9k*. 

Using the divergence theorem the following identity is easily established: 


(4.11) f r]{uhit, ■)) dx = [ r]{uh{t+, ■)) dx 

jR<i JRd' 


+ J2 ^-FMdx- 


kGTX 


«eT/* 


I Ok' 


• irds. 


We consider the entropy rj = for even p > 4. If we re-write the surface integral 
over 9 k* to include contributions from both sides of the interface, the last term on 
the right hand side is equal to 


1 


1 


.Me) 


/ iF{u[^’‘>)-F{u[^>))-i^ds = - 


- «er," '^9., 

Using the pointwise estimate 


kGTA “)> 


F'{f)-vd(,ds. 




\i=0 


\i^0 


we can bound 


P'^h ' P'^h ‘ C' C 'I 

(4.12) F\f)-ndf<C |^'(^)|d^<_|(4«<=))P + (4'^))i>|. 
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Consequently, re-distributing the surface integrals to include contributions only from 
the interior face, one has 


(4.13) 


1 

2 


i 

.er," 


< 



where a trace inverse inequality has been used. 

We now consider the divergence term in (14.1111 . To that end, we first note the 
pointwise estimate 




1,1 7/1, 


whence 


f y-F{uh)dx<—Y / I^K)! da; 




(4.14) 


C 




KGTy 


Substituting (14.131) and (14.141) into (14.111) , we see that for some positive constants Ci, 
C 2 , we have 


f Uh{t,xYdx<^ f f Uh{T,x)P dxdr + C 2 - 
jRa h Jtn aRrf 

Since we already have discrete stability ||uh(t“, •)||j;,p(Rd) < 00 for for n = 
0,1,..., A^, a simple Gronwall argument allows us to conclude 

sup I|u?,(t, •)||LP(Ra) < c, 

for all h and p satisfying (gH). Using an inverse estimate, the rest of the proof can 
be completed as in [nmu]. □ 

5 . Entropy Consistency. We must prove the weak consistency of scheme (12.31) 
with all entropy inequalities. This is formally stated in the following theorem: 

Theorem 5.1. Let Uh be the solution produced by scheme (12.31) . Then, for all 
convex entropies, and associated space-time entropy fluxes dSH), there holds 

(5.1) liminfM F(m^) • V^dx i > 0, V((> € C(i”((0, T) x R+). 

^^0 l^J(0,T)xRa J 

In the remainder of this section we prove Theorem l5.ll It is well known that it is 
sufficient to show dSI]) for all Kruzkov entropy pairs, i.e. convex entropies rjk = \u — k\ 
with fc £ R, and associated entropy fluxes m HI- Following mm we define a 
regularized Kruzkov entropy function 7]k,s, converging to r]k as i5 —>-0, uniformly on 
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compact sets, while maintaining convexity and uniformly bounded first derivatives. 
This allows us to define an associated regularized entropy flux which converges 
uniformly on compact sets to the Kruzkov entropy flux. Such a regularization can be 
accomplished with standard mollifiers. A more detailed exposition can be found, for 
example, in lainiiii]. For the remainder of this section we simply identify 7] = 7]k,5 
and F = for some fixed ^ > 0, which allows us to assume differentiability and 
that all derivatives 77 *-'^^ are bounded. After establishing (15.1|) for the regularized 
entropy pairs, we let <5 —>■ 0, and conclude theorem 15. ll for the Kruzkov entropy pairs 
by applying dominated convergence. 

Before we proceed with the proof, we give a few preparatory lemmas. First, note 
that in (HU) there is no term controlling the L^-norm of the flux divergence. Such 
terms arise in a straightforward manner from the streamline-diffusion stabilization [H 
0 [TH] . The next lemma shows that we can control the flux divergence even in the 
absence of the streamline-diffusion stabilization. We shall need this result to complete 
our proofs. 

Lemma 5.1. Let Uh he the solution produced by scheme (lOl) . Then there exists 
a constant C > 0, independent of h and Uh, such that 


(5.2) 



h^^\V ■ f{uh)\^ dx < Ch^ 


where j3 is defined in IST71 . 

Proof. Split the integral in (15.21) into integrals over := : |Vm?i(x)| < h~3 } 

and := Kn{x : |Vu?i(a;)| > }. Using (12.9|) we have 


h^^y^ f \V-f{uh)\'^dx<Ch^^yf \Vuh\^dx<Ch^, 

„ Jk< „ Jk.< 


and 


h^^ E/ IV ■ f{uh)\‘^dx<ch'^^ [ \y ■ fiuh)\\yuh\dx 

<Ch'^^yf \V-f{uh)\\Vuh\‘^hidx 

n,K 

< C/ 72 / 3+1 f |v./(n,)|dx 

n,K ^ 

< C/i2/3-¥ ^ f e{uh)dx 

<Ch^y f £{uh)\Vuh\‘^ dx. 


Here we have used the definition of the artificial viscosity ( 1 ^ in the second-but-last 
inequality, and again the standard inverse estimate (14.61) in the last inequality. The 
lemma now follows from L2_g|;a]3jii^y (HU). □ 

Lemma 5.2. Let Uh be the solution produced by scheme (El, and let eiuh) he 
defined as in dSJl). Then there exists a constant C > 0 such that for small enough h 
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Proof. Using the definition (I2.7L we have 

f £{uh)dx = Y^ I h^\V ■ f{uh)\dx + Y^ j \ dx 

n,K ^ n,K n,K. 

- fe / ■ fiuh)\^ dx\ f dx^ 


J, 

2 


The lemma follows from shape-regularity (12.IL Lemma 15.11 and the L^-stability re¬ 
sults 611), 63D. □ 

Recall that we have assumed /3 > 5 , so the left-hand side of (15.31) tends to zero 
as h —0. We are now ready to prove the main result. 

Proof, (of Theorem l5.ll) We start by considering (13.21) for the regularized Kruzkov 
entropy and general non-negative 4> with compact support in (0, T) xK.'^. Furthermore, 
using the projection (13.11) . we obtain from (13.61) . noting (13.7L (13.81) . 


/ F{uh) ■ V(/)dx > Af^^{uh,TTh{v'(l^) - ??'<(>) +M^^{uh,v'(l))- 
Split the shock-capturing term as follows: 

f e\Vuh\'^r]"{uh)4>dx+ '^ f eVuh-V(j)r]'{uh)dx\ . 

\n,K n,K, / 

Clearly, the first term is non-negative, while the second term vanishes as h —> 0, i.e. 

/ eS/uh ■ V4>r]'{uh) dx < C | f elVit/ipdxj ( X] / 

vanishes by L^-stability 63) and Lemma [521 Now define Crj^h '■= ~ l'4‘^ 8 'iid 

write 

N°^{uh, e^^h) = X |y ^ ■ f('^h)en,h + J (^f - dsj . 


Using the definition of the shock-capturing operator 63, we finally obtain by the 
estimates provided in [T3] (see (3.8c) in [T31) 


n,K Jk ^ 

(5.4) < [ s{uh) {\Vuh\^ + 1) dx 


where (14.61) and the uniform L°°-stability. The remaining integrals are bound by in¬ 
stability and Lemma|n 21 and we can conclude that \Af^^{uh, e^,h)\ —>■ 0 , as h —>• 0 . 

We have thus proved Theorem 15.II for the regularized Kruzkov entropy pair rjs.k, 
Fs,k- Using the uniform convergence of rjs^k and Fs^k as i5 —>■ 0 on compact sets to rjk 
and ifc respectively, the entropy inequality 63 follows by dominated convergence. □ 
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6 . Conclusion. We have considered a class of space-time DG schemes for scalar 
hyperbolic conservation laws which do not use any streamline diffusion stabilization. 
We have shown that a shock-capturing term, somewhat simplified compared to similar 
schemes m, provides sufficient stabilization. For all orders of polynomial approxi¬ 
mation, the scheme admits a bound on the solution in L°“([0,T] x K.'^), and is weakly 
consistent with all entropy inequalities. 

Moreover, the scheme is consistent with the initial data. The proof is omitted 
here, as it clearly does not depend on streamline-diffusion stabilization. See m and 
the references cited therein for details on consistency with the initial data. 

Using DiPerna’s theorem [7] we can conclude that the scheme converges to the 
unique entropy solution of cu for all polynomial orders. Similar convergence proofs 
for nonlinear conservation laws have thus far only been available for schemes that use 
streamline diffusion stabilization. 
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Appendix A. Proof of Lemma l4.ll First, we review some properties and theo¬ 
rems in polynomial inequalities. More specifically, we introduce the Kneser inequality 
and its extension to multivariate polynomials which is the main interest here. Then 
some properties in the equivalency of norms in finite dimensional spaces are revisited. 
Then we are ready to give a short proof of Lemma 14.11 
Recall the classic Kneser inequality [2j p. 260]: 

Theorem A.l. Suppose f = gh where g G P™ and h G where P™ 

denotes the polynomials of order m of one complex variable and complex coefficients. 
Then 

(A-l) Il5ll[-l,l]ll^ll[-l,l] ^ 2^n,mCn,n-m\\f\\[-l,l], 

where by || • ||[-i,i] we mean the uniform norm on the interval [—1,1] of real numbers 
for continuous functions, and the coefficient Cn,m is defined as 

Cn,m ■■= 2™ n (1 + COS . 

One might simplify the complicated form of the coefficients Cn,m with this corollary 
from [3J p. 261]: 

Corollary A.l. Under the assumptions of Theorem \A.l[ the following holds: 


(A.2) 


/ ||g||[-ij]||fe||[-i,i] N 1/" < 3 3 
^ II/II[-L1] ^ 


By a change of variables, inequality (IA.1|) can be easily extended to any real 
interval J = [a, 6], i.e. 

(A.3) \\g\\j\\h\\j<KMj- 

Let us call Kn Kneser’s constant, which only depends on n. Using (IA.2p . we can 
set Kri = 3.3” . Also note that by the definition of the uniform norm we may write 
\\-\\j = \\-\\L^ijy 
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Here we are interested in the case where the polynomial coefficients are real, and 
we seek the extension of Theorem lA.ll to general multivariate real polynomials of total 
degree n in a convex domain J C We can prove the following theorem, based on 
o proof of Theorem 2.2]: 

Theorem A.2. Assume J C is convex and f is a polynomial of d variables of 
total degree n defined on J such that f = gh where the factors g and h are polynomials 
of total degree less than n. Then the following holds 

(A.4) ||5|U”(J)II^I|l“(J) < Kn\\f\\L°°{J), 

where Kn is the same as the constants defined in (IA.1|) or (IA.2I) . 

Proof Let us assume the xq^xi S J are the maximum points of g and h respec¬ 
tively, i.e. 

llffl|L“(j) = l5(a:o)|, Whh’^ij) = \h{xi)\. 

Define the line segment ^ := {x = txi + {t — l)xo; t G [0,1]}. Due to convexity of J, 
one knows that S G J. So the Kneser inequality holds on this line 

llffl|L“(S)ll^l|L“(S) < ^n||/||L~(S)- 
The special property of segment S, xq,xi € S, gives 

On the other hand it is clear that ||/||l~(S) < II/I|l“(j), so ( |X4t holds. □ 

Also we will use the following equivalency of norms for a general function / 
in a finite (iV-) dimensional space Xat and 1 < 9 < p < 00 , 

A-5) |j|./]-./, ll/IU-w < ll/IU-(J) < 

where the left inequality comes directly from Holder’s inequality and the right in¬ 
equality can be found in [3] p. 140]. Note that by the dependency on X^r we mean 
the dependency on the type of this space like d and N (or, in our case, equivalently 
the polynomial total degree n). 

Now we are ready to present the proof of Lemma 14.11 

Proof. [Proof of Lemma 14.1| Let us consider J = k in the result showed before. 
By definition of the P[^ projection (13.11) one has 


f Vv ■V{TTh{v^ ^))da; 


Vu • V(u^ ^)dx=(p—1) 


^ dx. 


^ K J K ^ K, 

The goal of the proof is to find a lower bound for the |||Vupu^’“^||ii(K) for even integer 
p as the following 

(A.6) II|Vu|V-2|Ui(„) > C|||Vul2]Ui(,)]]u^’-2]U„.(,) 

Here ]Vz;p and are polynomials; since v G we know that 


Vv € 


yP-^ G G 
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Using the equivalence of norms {L^ and L°°) as in (IA.5I) gives 

ll|v..|V-"||yM > c(pp!-2) IIF"IV-"|U.,., 

Next we apply the Kneser inequality in the form of (IA.4I) and again norm equivalency 
(IA.5P to get 


lllVul 


2„P-2| 




1^1 


Kpg_2C{rp<i-^) 


> 


1 


Kpg_2C{VP<l-^) 


iiivunuo.(.)iiu^-2|u=o(.) 

iiivunuq«)iiu^-2|u=.(«). 


So the lower bound constant in (|A.6I) is 


Kpg_2C{VP‘i-^) 


; which yields 


- - [ VvV{T:h{vP i))da;. 

P ^ J K 


□ 

For brevity we show the constant 
as C{q,p). 


Kpg_2C{VP^-^) 

p-l 


by its two effective parameters 
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